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Predictions for the transverse single spin asymmetry (SSA), An, are given for the inclusive pro- 
cesses ip' —¥ h X and lp^ — > jet + X, which could be measured in operating or future experiments. 
These estimates are based on the Sivers distributions and the Collins fragmentation functions which 
fit the azimuthal asymmetries measured in semi-inclusive deep inelastic scattering (SIDIS) processes 
(£p^ — > £' hX). The factorization in terms of transverse momentum dependent distribution and 
fragmentation functions (TMD factorization) - which supplies the theoretical framework in which 
SIDIS azimuthal asymmetries are analyzed - is assumed to hold also for the £p — > hX inclusive pro- 
cess at large Pt- A measurement of A m would then provide a direct test of the validity of the TMD 
factorization in this case and would have important consequences for the study and understanding 
of SSAs in pp' — ¥ h X processes. 

PACS numbers: 13.88. +e, 13.60.-r, 13.85.Ni 



I. INTRODUCTION 



Transverse single spin asymmetries (SSAs) in semi-inclusive deep inelastic scattering (SIDIS), IN — > i'hX, have 
been measured by HERMES and COMPASS @-[|. A lar ge amount of data is still being analyzed by these 

Collaborations and new results are expected soon from the JLab experiments at 6 GeV. A rich program focused 
on azimuthal asymmetries, as a way of probing the internal nucleon structure, is planned for JLab operating at an 
upgraded energy of 12 GeV and for the future electron-ion (EIC) or electron-nucleon (ENC) colliders, which are under 
active consideration within the hadron physics scientific community (see e.g. Ref. 0] for a short up-to-date overview). 

These SIDIS SSAs are interpreted and discussed in terms of unintegrated, transverse momentum dependent, dis- 
tribution and fragmentation functions (shortly, TMDs) . In particular the Sivers distributions [l(| [HI and the Collins 
fragmentation functions [l2| have been extracted fl3l - flq from SIDIS data, and, thanks to complementary information 
from Belle on the Collins function [l9|, [2(| , a first extraction of the transversity distribution has been possible [U [22| ■ 

All these analyses have been performed in the 7* — p cm. frame, within a QCD factorization scheme, according to 
which the SIDIS cross section is written as a convolution of TMDs and elementary interactions: 

da ^e'hX = J- f q/p ( x> k±] Q 2 ) ® da e ^« ® D h/q (z,p ± ; Q 2 ) , (1) 

where k± and are, respectively, the transverse momentum of the quark in the proton and of the final hadron with 
respect to the fragmenting quark. At order k±/Q the observed transverse momentum, Pt, of the hadron is given by 

P T = k ± + zp ± . (2) 

There is a general consensus [23l - [27j that such a scheme holds in the kinematical region defined by 

P T ~ k± =s Aqcd < Q ■ (3) 

The presence of the two scales, small Pt and large Q, allows to identify the contribution from the unintegrated 
partonic distribution (Pt — k±), while remaining in the region of validity of the QCD parton model. At larger values 
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of Pt other mechanisms, like quark-gluon correlations and higher order pQCD contributions become important [27l — 
[29j | . A similar situation [24l [26l. [30l - [35| holds for Drell-Yan processes, AB — > l + t~X, where the two scales are the 
small transverse momentum, qt, and the large invariant mass, M, of the dilepton pair. 

The situation is not so clear for processes in which only one large scale is detected, like the inclusive production, at 
large Pt, of a single particle in hadronic interactions, AB — > CX. However, the most striking and large SSAs have 
been [36M39[ | and keep being measured (4pl - [i3 | in these cases. The TMD factorization for these processes was first 
suggested in Refs. Hoi. Illll a nd adopted in Refs. (45l - l47j to explain the large single spin asymmetries observed by the 
E704 Collaboration J37|, |39[ . The same approach led to successful predictions [48|, |49| for the values of An measured 
at RHIC [13]. 

Alternative approaches to explain the origin of SSAs, linking collinear partonic dynamics to higher-twist quark-gluon 
correlations, were originally proposed in Refs. [5l| - l55j and phenomenologically adopted in Refs. (56l-l59|. These two 
approaches, the TMD factorization and the higher- twist correlations, have been shown to be somewhat related [60l.l6l| 
and consistent with each other (3^, HH [62j . 

However, a definite proof of the validity of the TMD factorization for hadronic inclusive processes with one large scale 
only is still lacking. Due to this, the study of dijet production at large Pt in hadronic processes was proposed [63j-l66|, 
where the second small scale is the total qt of the two jets, which is of the order of the intrinsic partonic momentum 
k±. This approach leads to a modified TMD factorization approach, with the inclusion in the elementary processes 
of gauge link color factors [67t - l69j . 

In this paper we propose a phenomenological test of the validity of the TMD factorization in cases in which only 
one large scale is detected, by considering SSAs for the lp* — > hX process, with the detection, in the lepton-proton 
cm. frame, of a single large Pt final particle, typically a pion. The final lepton is not observed; notice, however, that a 
large value of Pt implies, at leading pcrturbative order, large values of Q 2 . Such a measurement is the exact analogue 
of the SSAs observed in the pp^ —> hX processes, the well known and large left-right asymmetries An [36|-[44j. We 
compute these SSAs assuming the TMD factorization and using the relevant TMDs (Sivers and Collins functions) as 
extracted from SIDIS data. 

Such a choice is natural for the Collins function, which is expected to be universal (70l . [7l| . The Sivers distribution, 
instead, is expected to be process dependent as it is originated by final (or initial, depending on the process considered) 
state interactions, which also model the gauge links necessary for its correct gauge invariant definition [2^, [72|, [73J. 
However, these final state interactions should be the same in usual SIDIS processes and in the process considered 
here. 

A similar idea of computing left-right asymmetries in SIDIS processes, although with different motivations and 



still demanding the observation of the final lepton, has been discussed in Ref. [74| . A first simplified study of An in 
ip^ — > h X processes was performed in Ref. [75j]. The process was also considered in Refs. [7(1 [77j in the framework 
of collinear factorization with twist-three correlation functions, obtaining anomalously large asymmetries with a sign 
opposite to that of the corresponding asymmetries in pp processes. 

The plan of the paper is the following: in Section [TT| we present the formalism for the study of SSAs in a TMD 
approach for both the p^£ — > hX and the p^i — > jet + X processes; in Section Hill we show our numerical estimates 
of the contributions of the Sivers and Collins effects to An, based on the present knowledge of TMDs, for several 
different kinematical setups and discuss their phenomenological aspects; finally, in Section HVl we give some comments 
and conclusions. Technical details on the full noncollinear kinematics are given in Appendix [Al while the calculation 
of the hclicity amplitudes is worked out in Appendix iBl The complete expression of An for the process p> (' — > h X. 
including all TMD contributions at leading twist, can be found in Appendix [Cl 

II. FORMALISM 

A. Large Pt hadron production 

We propose to study single spin asymmetries for the the process p^i — >• h X in close analogy to the study of the SSAs 
for the process p^p — > h X, assuming the validity of the TMD factorization. The cross section for this process can then 
be written as a particular case of the general treatment, in a factorized scheme, of the (^4, Sa) + (B, Sb) —> C + X 
large Pt inclusive polarized process [H, l7§ - l80j : 

Ehda ( P .s )+ e^ + x , dxd2 



d 3 P h 



= E / 1fi d ' T 2 dZ 2 d2fc ^ d3 P± S (PX ■ P' q ) J(P±) 6(S + i+u) 
J ID 7T X Z S 

<MA} 

x |< P a? \ Mk,\^\ ^V^Px) \ . 
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FIG. 1: Kinematical configuration and conventions for the p^£ — > hX process. 



which can be shortened, with obvious notations, as 

da; dz 



dcr- 



E 

9,W 



16 TT 2 X Z 2 S 



d 2 k ± d 3 P± 5{p ± ■ p') J(p ± ) 5(s + t + u) z, k ± , P± ) , 



(5) 



where is the term in curly brackets of Eq. (j4j. 

Let us recall the main features of these equations. 

• We consider the collision of a polarized proton (or, in general, a nucleon) in a pure transverse spin state S with 
an unpolarized lepton, in the proton-lepton center of mass frame. The proton p moves along the positive Z cm 
axis and hadron h is produced in the (XZ) cm plane. We define as transverse polarization for the proton the 
direction, often using the notation f and 4- respectively for protons polarized along or opposite to Y cm . The X cm 
axis is defined in such a way that a hadron h with (Ph)x cia > is produced to the left of the incoming proton. 
The transverse momentum is denoted as Pt- This kinematical configuration is shown in Fig. [TJ Results for the 
case of leptons moving along the positive Z cm axis (£p^ — > hX) will also be discussed in the paper. 

• The notation {A} implies a sum over all helicity indices, x and z are the usual light-cone momentum fractions, 
of partons in hadrons (x) and hadrons in partons (z). k± and p_|_ are respectively the transverse momentum of 
the parton q with respect to its parent nucleon p, and of hadron h with respect to its parent parton q. p' q is the 
three-momentum of the final fragmenting parton; it can be expressed in terms of the integration variables and 
the observed final hadron momentum. We consider all partons as massless, neglecting heavy quark contributions. 
Full details can be found in Ref. (80j and useful expressions are given in Appendix [X\ 

• With massless partons, the function J is given by [48{ 



J(P±) = 



E h 



■Pi 



4(P,2 



■pi) 



(6) 



In the kinematical regions which we shall consider J is close to 1. 
• p 9 / P y is the helicity density matrix of parton q inside the polarized proton p, with spin state S. f q / p ,s{ x i k±) 



is the distribution function of the unpolarized parton q inside the polarized proton p. 

q/p,s f 
P\ q ,\' q Jq/P: 

angle of k± 



P\ X fq/p,s( x ik±) are directly related to the leading-twist TMDs, with a dependence on 



The products 
the azimuthal 



The M 



A , A g \ A , A , 



's are the helicity amplitudes for the elementary process ql ~ > q£, normalized so that the 
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unpolarizcd cross section, for a collincar collision, is given by 

A„,A f 

At lowest perturbative order q£ — > qi is the only elementary interaction which contributes; notice that, in 
the presence of parton intrinsic motion, it is not a planar process in our chosen frame and depends on the 
intrinsic momenta, including their phases. Neglecting lepton and quark masses there are two independent 
hclicity amplitudes: 

M ++ . i++ (s,i,u,k±) = M*_.__ = -8ire q aj e ltpi = e iipi (8) 
M + _ ;+ _(s,t,u,fe_L) = M1+.-+ = STre q ale^ 2 = M 2 ° e lip2 , (9) 

where ifi^ are phases explicitly given in Appendix 151 Eqs. (|B8[) and (|B9|) . 

D\ h \I l {z,pj_) is the product of fragmentation amplitudes for the q h + X process 

q ' q 

n ^y ^/Va^Va^;- do) 

where the J stands for a spin sum and phase space integration over all undetected particles, considered 

J A', A x 

as a system X. The usual unpolarized fragmentation function D h / q {z) 1 i.e. the number density of hadrons h 
resulting from the fragmentation of an unpolarized parton q and carrying a light-cone momentum fraction z, is 
given by 



D h/q 



(*) = 5 E / dVi^j^frpJ- (11) 



We shall only consider the case of spinless final particles (X h = 0), in particular pions. In general D\ \i {z,p ± ) 
depends on the azimuthal angle of h around the direction of motion of the fragmenting polarized parton j80| . 



We compute the SSA: 



_ Ao-HP t ) - d^(Pr) _ Aa\P T ) - d^{-P T ) 
N dcrt(P T ) +do^(P T ) 2da un P(P T ) ' [ ' 



which can be measured either by looking at the production of hadrons at a fixed transverse momentum Pp, changing 
the incoming proton polarization from j" to J,, or keeping a fixed proton polarization and looking at the hadron 
production to the left and the right of the Z cm axis, see Fig.[TJ An is defined (and computed) for a proton polarization 
normal (N) to the production plane and a pure spin state (a pseudo- vector polarization St with \St\ = St = 1)- For 
a generic transverse polarization along an azimuthal direction <f>s (in our chosen reference frame) and a polarization 
St 7^= 1, one has: 

A{(j> s ,S T ) = S T -(px Pt)A n = S T sin(psA N . (13) 

Notice that if, according to the usual procedure in SIDIS experiments, one defines 

.sinfe _ 2 / d(p s [dcr((/> s ) - da{cj) S + ir)} sin^g 
TU ~ S T f 64> s [d*(<i>s) + M<t>s+*)] ' 



one simply has 



A™** = A N . (15) 



In order to compute An, Eq. (|T^|) . we need to compute [S(t) — £(!)] and [E(t) + £(!)], which can be done by 
performing the helicity sum in Eqs. (jH) and As our process is a simple particular case of (A, Sa) + (B, Sb) —> C+X, 
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the result agrees with Eqs. (82) and (86) of Ref. [8(|, simplified to the case in which particle B is a point- like lepton and 
the elementary interaction has only two independent amplitudes. Notice that several TMDs appear in the expression 
for j4jv; however, numerical evaluations show that the contribution of the Sivers effect is the dominant one. A modest 
contribution is given by the Collins function (coupled to the transversity distribution), while another contribution 
involving h\ T (sec Appendix [C]) is totally negligible. Considering only the Sivers and Collins effects, one has: 



E / t/^^ ^ d "P± S (p± ■ P' q ) J (P±) S(a + i + u) [E(t) - m] q ^ qi 

E / J X 2 A \ d2fe ^ d3 p± s (p±- ■ p q ) j (p^ s ( § + 1 + ") p(t) + m] q ^ qt 



with 



{A} 



[E(t)-S(4-)]^ = lA% /pt ( Xl k ± )co^ [|M?| 2 + |M°| 2 1 D h/q (z, P± ) 



+ h lq (x, k ± )M°M° A N D h/qt (z,p±) cos(0' + <f> h q ) (17) 
and (dropping negligible contributions from other TMDs j80| ) 

^[S(t) + S(4)]^=/ 3/p (x,^) [|M°| 2 + |M 2 T] D h/q (z, Px ). (18) 

{A} 

• The first term on the r.h.s. of Eq. (fTTj) shows the contribution to of the Sivers function A N f q / p t(x, fcjj [T3, 

El El, 

&fq/p,s( x > k± ) = fq/p,s(x, k±) - f q /p-s(x, k±) = A N f q/pt (x, k±) S T -(px k±) (19) 

= -2^jf^(x,k ± )S T -(pxk ± ), 

coupled to the unpolarized elementary interaction (<x i ( | AiTj 1 1 2 + \M% | 2 )) and the unpolarized fragmentation 
function D h / q (z,p±); the cos <f> factor arises from the St • (p x fej_) factor, the spin-transverse motion correlation 
of the Sivers function in the case of a normal spin direction with St = 1- 

• The second term on the r.h.s. of Eq. (|17l) shows the contribution to A n o f the unintegrated transversity distri- 
bution hi q (x,kj_) coupled to the Collins function A N D h / q t(z,p±) [l2ll8ll|. 

AD h/qr (z,p ± ) = D h/q t (z,p ± ) - D h/q i (z,p ± ) = A N D h/gf (z,p±) s q ■ (p q x p ± ) (20) 

= -^Hi q (z,p±) s q ■ (p x p ± ) , 
zmh H 

and to the transverse spin transfer elementary interaction (dcr"^ — dcr>+ oc M^M®). The factor cos((f>' +4> q ) arises 
from phases in the fcx-dependent transversity distribution, the Collins function and the elementary polarized 
interaction. cf>' is the azimuthal angle of the fragmenting quark (with 3- momentum p' q ) and <fi q is the azimuthal 
angle of p ± around the p' q direction (80j . Their expressions in terms of integration and overall variables can be 
found in Appendix |A1 

• The elementary interaction amplitudes are explicitly given in Eqs. ([5J and ([9]). Notice that the elementary 
Mandelstam variables s,t,u are computed taking into account the full kinematics, and thus depend on the 
transverse momenta. 

• A final issue which needs to be clarified concerns perturbative QCD corrections. Our proposed process involves 
TMDs coupled to lowest order perturbative interactions and is driven by a large angle elementary electromagnetic 
scattering, ql q£. Some QCD effects, like soft gluon emissions, are taken into account in the TMDs, as the 
emission of soft gluons builds up intrinsic partonic motion. Higher order pQCD corrections due to genuine hard 
QCD processes, like ql — > qig orgf-> qql are not included in our computation of An- These contribute at 
order a s to the cross section and can be neglected at large Q 2 values; moreover, one should notice that events 
induced by these hard pQCD elementary interactions result in final states with two fragmenting partons, i.e. two 
jets, and could be experimentally excluded. However, these pQCD corrections might be of some relevance and 
difficult to disentangle at HERMES, COMPASS or JLab energies. 
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B. Large Pt jet production 

We consider also the most interesting case of SSAs for the inclusive process p^ £ — > jet + X. Although it is a 
difficult process to detect experimentally and might require future higher energy and luminosity machines, it would 
certainly give the most direct access to the Sivcrs effect, as the lack of any fragmentation mechanism forbids other 
contributions. Even more difficult, the observation of both a jet and a final hadron inside the jet (with a measurement 
of its transverse momentum p ± ) , would allow a direct detection of the Collins effect [82| ■ 

In the case of the p^ £ — > jet + X process, with no observation of a single final particle, Eq. (j4j) simplifies to: 



d3R 



(l.C 



16 ir 2 x s 



d 2 k±_ 5(s +i+u) 



E 

q.{\\ 



(21) 



while Eq. ([TBI becomes: 



A) ct 
N 



E 

9,{A} 



E 

?,{A} 



d.r 



16 TT 2 X S 



d 2 k ± S(s + t + u) [S(t)-S(|)]| 



ql^tql 
ct 



dx 



16 TT 2 X S 



d 2 k ± s(s +i+u) [s(t) + m)]£~* 9t 



J jot 



(22) 



In this case the kinematics is very simple and is shown explicitly in appendix I A 21 For a generic azimuthal direction 
bs of the transverse spin St, the Sivers function, Eq. (|19|) . can be written as: 



A /,/pt (x,k±) S T ■ (px fcx) = A / g/pt (x, fc±) ^sm</> s — - cos0 s ^ 

= & N fq/ P i (x,k±) sin(0 s - 0) , 



k x . 



(23) 



and the £ kernels in Eq. 



are 



£ [S(t) - S(|)]^ 9£ = 5 A% pt (x, fc ± ) sin(0 5 - < 

{A} 

£[£(t) + ^)]?^ = /,/*(*, fc-0 [iMrT + l^l 2 

{A} 



|M°| 2 



(24) 
(25) 



The elementary amplitudes are the same as given in Eqs. © and ©• 



III. ESTIMATES FOR ^jv 



We have computed the SSA, An, as defined in Eq. (|12[) or (Ti"4"|) . for the large Pt production of pions and jets in 
p'£ — > h X and — > jet + A processes, according to the expressions given, respectively, in Eqs. (|T6l) - (fT8)) and in 
Eqs. (JUJ, (J2U) (with c/> s = tt/2), and ([23]). 

Analogous results for the case of leptons moving along the Z cm axis, £p^ — > h (jet) +X, in the same chosen hadronic 
frame (that is, keeping fixed the definitions of xf = 2Pl/\/s and of the f, J, transverse polarization directions) can 
be easily obtained using rotational invariance: 

Af^ et)+x (x F ,P T ) = -A^ h ^ +x (-x F ,Pr). (26) 

We have used the Sivcrs distributions as parameterized and extracted - from SIDIS data - in Rcf. fl4|: even if 
the Sivers functions, being related to final state interactions [72j], are expected to be process dependent |23|. they 
should be the same in SIDIS and the (related) processes considered here, which all originate from the same q£ —> q£ 
elementary interaction and subsequent quark fragmentation. Similarly, we have used the transversity distributions 
and Collins functions as parameterized and extracted in Ref. [22j |. The un pola rized parton distribution functions 
(PDFs) and fragmentation functions (FFs) are taken respectively from Refs. j83[ and [84j . 
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FIG. 2: Estimates of A N vs. x F for the p 1 1 -> n X process at HERMES (y/s ~ 7 GeV). Left panel: Sivers effect at P T = 1.5 GeV; 
central panel: Sivers effect at Pt = 2.5 GeV; right panel: Collins effect at Pt = 2.5 GeV. The computation has been performed 
according to Eqs. (116[) . (|17[) and (JTHJ) of the text, adopting the Sivers functions of Ref. [3] and the transversity and Collins 
functions of Ref. [22t|. as extracted from SIDIS and e + e~ data, the unpolarized PDFs of Ref. and the FFs of Ref. In 
the left panel we also show, for charged pions, the statistical uncertainty bands coming from the extracted Sivers functions 



Our results are given for the kinematical configurations of HERMES, COMPASS, JLab at 12 GeV, and a hypo- 
thetical ENC future machine operating at an energy yfs = 50 GeV. For hadron production, the Sivers and Collins 
contributions are shown separately. We plot Am as a function of Xf at fixed Pp values; these should be chosen as the 
hard scale of the process, ensuring a large momentum transfer in the hard scattering, say Q 2 > 1 GeV 2 . In collinear 
cases, at LO, it might suffice to have Pt > 1 GeV; however, with TMD factorization, one has to be more careful, as 
Pt might be partially generated by intrinsic k±. We have checked that a value of Pt = 2.5 GeV corresponds to a 
safe Q 2 > 1 GeV 2 region in the whole range of xf, while Pt = 1.5 GeV implies a safe Q 2 region only for backward 
production, xf ^ 0. We give predictions for these two values of Pp. 

Notice also that for positive xf the minimum of x is given, roughly, by xf- This implies that for xf > 0.2 - 
0.3 we should employ the parameterizations of the Sivers and transversity functions in a region where they are not 
constrained by SIDIS data. For this reason we will give our theoretical estimates of Am only up to xp — 0.2. On the 
other hand, for negative xp the minimum of x is controlled by the ratio xp = 2Pt/\/~s, implying that at moderate 
cm. energies (i.e. yfs ~ 10 - 20 GeV) and Pt — 1 - 2 GeV, we are sensitive to the valence region of the polarized 
proton, i. e. the region where the Sivers (and the transversity) functions reach their maxima. 

Let us comment in details our results. 

• We first stress some aspects peculiar to the p^£ -> ftl process. As in SIDIS processes at leading order accuracy, 
only one partonic subprocess, q£ — > qi, is active, with a simple 1/t 2 dependence (a much simpler dynamics 
than in the pp — > hX case). However, since the lepton plane is not identified (we do not require the detection 
of the outgoing lepton), one cannot access, separately, the Sivers and the Collins effects. Nevertheless, in the 
backward region (w.r.t. the proton direction) the variable |u| becomes smaller and so does the partonic spin 
transfer cross section cx M^M® [see Eqs. ([H]) and (JSJ)], entering the Collins contribution to Am [second term 
on the r.h.s. of Eq. (|17[) ]. This implies a strong dynamical suppression of the Collins effect (reinforced by the 
integration over the azimuthal phases) at largely and moderately negative values of xf , leaving active mainly 
the Sivers contribution. Notice that, contrary to what happens in the pp -> ftl process, no u-channcl in the 
partonic process is present; moreover the variable t strongly depends on <p, the azimuthal phase of the Sivers 
effect [first term on the r.h.s. of Eq. (fTT)) ]. 

• In Fig. [2] we present our estimates, separately, for the Sivers and Collins contributions to Am at HERMES 
kinematics. More precisely, we show the Sivers effect at Pt = 1.5 GeV (left panel) and at Pt = 2.5 GeV 
(central panel) and the Collins effect at Pt = 2.5 GeV (right panel). The Collins effect at Pt = 1-5 GeV (not 
shown) is almost negligible in the kinematical region considered. For charged pion production at Pt =1.5 GeV 
(left panel) the statistical uncertainty bands as resulting from our fit [bi| are also shown. 

The largest Am values obtained correspond to the x region (of the polarized proton distributions) where the 
Sivers functions, for u and d quarks, reach their maxima. It is interesting to note that the sizable value of Am 
for 7T _ production (larger than the corresponding Sivers contribution to Ajjp in SIDIS) is due to the dominance 
of the d quark with a small contamination from the u quark. This is related to the fact that the light-cone 
momentum fraction z is always bigger than the maximum between |a;j^ | and xp, implying, at moderate and 
large \xf\, a dominance of the leading fragmentation functions. 
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FIG. 3: Same as in Fig.[2]but for COMPASS kinematics (y/s ~ 17 GeV). 



In Fig. [3] we show the analogous results for COMPASS kinematics. Again at Pp = 1.5 GeV only the Sivers 
effect gives a sizable contribution (left panel), while the Collins effect (not shown) is compatible with zero. At 
Pt = 2.5 GeV the Sivers effect (central panel) dominates only in the backward region, while in the forward 
region the Collins effect (right panel) becomes sizable. For charged pion production at Pt =1.5 GeV (left panel) 
the statistical uncertainty bands as resulting from our fit flij are also shown. The main difference w.r.t. An 
for HERMES kinematics at P T = 1.5 GeV (compare Figs. [Hand H left panels) is that at the larger COMPASS 
energy (y/s ~ 17 GeV) the valence region for the polarized proton, where the Sivers functions reach their 
maxima, starts dominating at larger xp. 

In Fig. U] we show our results for ENC kinematics at y/s = 50 GeV. For Pt = 1.5 GeV (left panel) only 
at the upper range of the safe xp values (i.e. xp < 0) the Sivers effect gives a sizable contribution, of the 
order of few percent (the Collins effect is once again negligible). At Pp = 2.5 GeV both the Sivers (central 
panel) and the Collins (right panel) contributions are comparable and sizable around xp ~ 0.2, therefore hardly 
distinguishable. This can be understood because at such Pp and energy values the valence region for the 
polarized proton dominates only for xp > 0, where both effects are active (see our comment on the suppression 
of the Collins effect for negative if at the beginning of this Section). 
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FIG. 4: Same as in Fig. [2] but for ENC kinematics at y/s = 50 GeV. 



• In Fig. [5] we show analogous estimates of the Sivers contribution to An for JLab kinematics at the upgraded 
energy £/Lab = 12 GeV, corresponding to a cm. energy y/s ~ 4.9 GeV. Again, in order to guarantee a sufficiently 
large momentum transfer we show results at Pt = 1.5 GeV vs. xp < 0.1. The results are comparable to the 
corresponding estimates for HERMES kinematics, see Fig. [3] (left panel), with large asymmetries (in size) for all 
pions. Given the lower cm. energy, however, cross sections are in general smaller than those for HERMES and 
COMPASS kinematics. Larger values of Pt are probably out of reach at Jlab, while the Collins contribution is 
again negligible in the xp region considered. 

• Finally, in Fig. [5] we show estimates of the Sivers contribution to An for the process p^£ — > jet + X for ENC 
kinematics (y/s = 50 GeV) at Pp = 1.5 GeV (left panel) and Pt = 2.5 GeV (right panel) vs. xp. The results 
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FIG. 5: Same as in Fig. [2] left panel, but for JLab kinematics at £"Lab = 12 GeV (\/s ~ 4.9 GeV). 




FIG. 6: Estimates of An vs. if for the £ — > jet + X process and for ENC kinematics at y/s = 50 GeV. Left panel: Sivers effect 
at Pt = 1.5 GeV; right panel: Sivers effect at Pt = 2.5 GeV. The computation has been performed according to Eqs. (|22l) . 
(|24[) and_ j25[l of the text, adopting the Sivers functions of Ref. [3], as extracted from SIDIS data, and the unpolarized PDFs 
of Ref. "~ 



are similar, both in size and shape, to the corresponding ones for neutral and positive pions, see Fig. 21 left 
and central panels (notice the different scale). The asymmetry is almost negligible at negative xf and becomes 
sizable only at the upper range of the safe Xf values. We have found that Am becomes even smaller at larger 
cm. energies. 



IV. COMMENTS AND CONCLUSIONS 



In this paper we have presented a phenomcnological study, based on the assumption of TMD factorization, of 
transverse single spin asymmetries for the inclusive production of large Pt pions and jets in lepton-proton collisions, 

I — > h (jci)+X. These asymmetries, measured in the lepton-proton cm. frame (since the final lepton is not observed, 
the 7*-proton cm. frame cannot be reconstructed), should involve the same TMD distribution and fragmentation 
functions which contribute to the transverse azimuthal asymmetries measured by the HERMES and COMPASS 
collaborations in the last years in semi-inclusive deep inelastic scattering. 

Using best-fit parameterizations of the TMD functions extracted from HERMES, COMPASS and Belle data, we have 
shown that, in the kinematical regions where our perturbative approach should be reliable, the asymmetries dominantly 
arise from the Sivers effect in the distribution sector and marginally from the Collins effect in the fragmentation 
sector (not present in the case of jet production). We have presented results for several kinematical configurations 
corresponding to present experimental setups (HERMES and COMPASS), to the forthcoming 12 GeV upgraded 
JLab setup and to a class of lepton-proton (ion) colliders (ENC) currently under active study in the QCD and hadron 
physics community. These results show that for pion production the Sivers Am can be sizable, at least for HERMES, 
COMPASS and JLab at 12 GeV kinematics. For pion and jet production at typical energies of the proposed ENC 
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colliders the asymmetries are much smaller and become larger only at the boundary of the safe kinematical regions, 
where, for pions, both the Sivers and the Collins contributions play a role and the two mechanisms cannot be 
disentangled. 

The measurement of these predicted asymmetries allows a test of the validity of the TMD factorization, largely 
accepted for SIDIS processes with two scales (small Pt and large Q), but still much debated for processes with only 
one large scale (Pt), like the one we are considering here. A test of TMD factorization in such processes is of great 
importance for a consistent understanding of the large SSAs measured in the single inclusive production of large Pt 
hadrons in proton-proton collisions. 

We stress once more that our predictions refer to large Pt production, in the lepton-proton cm. frame, at leading 
perturbative order. It implies that, in order to compare experimental data with our results, one has to select large 
Pt, single-jet events, excluding those events containing a second jet in the opposite hemisphere w.r.t. to the primary 
observed jet (containing the final observed hadron). This should avoid large Pt jets (or hadrons) coming from next-to- 
leading order partonic processes (hard pQCD corrections). Although these requirements might correspond to smaller 
cross sections and difficult selection procedures, we believe that the relevance of testing TMD factorization in this 
simple process justifies efforts in this direction and motivates our work. 
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Appendix A: Kinematics 
1. Hadron production 



We work in the proton-lepton center of mass frame, with the incoming proton and lepton moving along the Z c 
axis and the outgoing hadron emitted in the (XZ) cm plane: 



p= ^(1,0,0,1) 

£=^(1,0,0,-1) 
Ph = {E h ,P T ,0,P L ) 



El 



Pi 



Pi, 



where s is the proton-lepton cm. square energy and where we have assumed all particles to be massless. 
kinematical variables for the elementary underlying process result in (kj_ = \kj_\, p± = \pj_\) 



Pq = 

t = 

Pa = 



2Xy/s 



, fcj 



Xy/s _ k\ \ 

2 Ixsfs) 



2 



(1,0,0,-1) 



Pq 



2z 

-P'a 



Pl 



(p t -p1,-p^p l -p1) 



(Al) 

(A2) 
(A3) 
The 

(A4) 
(A5) 

(A6) 
(A7) 



with being the intrinsic transverse momentum of parton q inside the parent proton and p^ being the intrinsic 
transverse momentum of the detected final hadron h with respect to the fragmenting parton q' . The expression for 
p' q has been obtained by requiring z to be the light-cone momentum fraction of the emitted hadron, z = P^ /p' q + 
as defined in the helicity frame of the fragmenting quark q' , which we will denote as S. With this kinematics, the 
partonic Mandelstam invariants are 



xs 
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t = 



2: 



\/ei-pi-(^^-) 



2k x ± (P T -p x ± ) - 2k v ± p v ± 



u = -— 1 
2z \ 



E 2 h -Pl + P L -P] 



(A8) 



Notice that the orthogonality between p q and p ± , explicitly guaranteed through the delta function S(p ± ■ p' q ) in 
Eq. (|16[) , allows us to fix one component of the vector p ± in terms of all the others; in particular it gives 



IPj 



= P T pl+P LP ] 



p l i = ±JPrPl + PlP1- {Pi? - (Pi) 2 



(A9) 



Similarly, the other delta function in Eq. (|16[) . S(s + t + u), can be used to perform the integration over the light-cone 
fraction z fixing 



E h 



2x^/s 



1 + x- 



kl 

xs 



Pj_ + I 1 - X 



2k x _(P T -p x ± ) - 2k v ± p v 



(A10) 



with p V j_ given by Eq. (|A9[) . 

The angle c/> q , which identifies the direction of p± around p' q , can be expressed in terms of the p^ components, p x ^_, 
p y ± and p\ , simply by noticing that in the helicity frame of parton q' (where the Z axis coincides with the direction 
of p' ) this angle is the azimuth of p ± , that is: 



Y 



(All) 



The helicity frame S of parton q' can be reached by performing two rotations, as explained in Appendix C of Ref. j80{ . 
in the following way 



Z = P Q 



1 



Y — Z cm x p„t — 



X = Y x Z = 



-.(P T -pl,-pi,P L -pi) 

(pI,Pt-p x ± ,o) 



{Pt-pI){Pl -pD, -pI(Pl-pI), -(p1) 2 -(Pt-p1) 2 ] 

VE 2 h -pl-(PL-pi) 2 ^W^pl 
where p qT is given by the transverse components of p' q in the center of mass reference frame, S: 



PqT 



VK-p1-(Pl-p1} 



? (p T -pl- p y ± ,o), 



(A12) 
(A13) 
(A14) 

(A15) 



and p v ± is fixed by the orthogonality condition of Eq. (|A9|) . 
By replacing Eqs. (|A12|) - (|A14|) into Eq. (|A11|) we find 



. . h pi P T 
sin <p„ = — — , 

q p± VeI-pI-(Pl-pI) 2 



P'a 



PI- V E 'k-Pl-(PL-Pl) 2 ' 

Alternatively, in terms of angles instead of components we can write 



(A16) 



cos 



Pt . 
sin < 

P± 

Pi i 

pj_ sin 6 



cos 8 



smt 



(A17) 
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where <j)',0' are the azimuthal, polar angles of p' q and 9± is the polar angle of p ± in our cm. reference frame. 
Finally, the cos(</>' + (f>g) azimuthal dependence of the Collins effect, sec Eq. (fT7|) . can be expressed as 



cos(^ + $) 
or, more simply, in terms of angles 



Pi- [El-pl-{P L ~piY] 



cos(</>' + = — sin 2 </>' - cosO^^- 
q p± sin ( 



(A18) 



(A19) 



2. Jet production 

The 4-momenta involved, in our reference frame and neglecting all masses, are 



P 



Pq 



(1,0,0,1) 



(1,0,0,-1) 



X\ s 



k I X \/~S k I 

k± , — 



2xy/s' 



p' q =P 3 = (Ej, Pt, 0, Pl,) 
so that the partonic Mandclstam invariants are given by 

s = xs 

t = 2P T k\ - x^ 
u = -^(Ej+P L ). 



2xyf S 



E* = P*+Pl, 



{Eo+Pl) 



(A20) 

(A21) 
(A22) 

(A23) 
(A24) 
(A25) 



Notice that there is no linear k v ^_ dependence in these variables and, as a consequence, in the elementary amplitudes 
Mi 2 - The delta function ensuring s + t + u = can be used to perform the integration on k^_ in Eq. (|22|) : 



Pt 



± 



V~s 1 (k y ± f 



Ej +P L V Ej +P L x 



(A26) 



which implies 



Ej+Pl 
2yfs 



4(fc^) 2 



yfs{Ej+P L ) 



(A27) 



Note that the term proportional to cos <ps m Eq. ([23f , being odd in k v ± , vanishes when integrating over k v ± , resulting, 
as it should, in a SSA proportional to sin 0s. 



Appendix B: Spinors and helicity amplitudes 



We compute the helicity amplitudes for the non- planar q£ —> ql process exploiting the well known spinor helicity 
technique (see, for example, Refs. (85l. l86l]). To be precise, we adopt the phase convention and gamma matrix repre- 
sentation of Ref. [86|; that is, our helicity spinors for a massless Dirac particle with 4- momentum k = (k°, k x ,k v ,k z ) 
and helicity ±1/2 are given by: 



u+(k) = v-(k) = 



/Vk+e-^/ 2 ' 
y/k = e l<t '' 2 



( 



\ 



U-(k) = w+(fc) 



(Bl) 
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where 



k x ± ik y 



k x ± ik v 



yj{k x ) 2 + (kv) 2 ' Vk+k~ 



fc± = k°±k z 



(B2) 



The two independent helicity amplitudes -Ma 3 ,a 4 ;Ai,a 2 for the q(ki, Ai) + £(k2, A2) — > 9(^3, A3) + £(fc4, A4) elementary 
lowest order QED interaction are given by: 



M. 



+ + : + + 



M*__._ 



= M*_ + ._ + = +2- 



[3 4] (12) 
[2 3] (14) , 



with 



Eq. (|B5P can be rewritten as [Hj]: 



V3) 



where = (fej + k./) 2 



COS( 



-i(<Pi+<f> j )/2 



2/c, • kj, and 



sin ( 



\J\ s i] \ fctkf 

With our kinematical configuration {<p2 = 7r) we obtain: 



Sij I fc^K^ 1 " 



0ij = 0ji + 7T 



(B3) 
(B4) 

(B5) 
(B6) 

(B7) 



M++.++ = M*__.__ = -8 7re g a * e-^ 31 e i(<?i3+ * 4 - 01+7r)/2 



M 4 



= M* 



q t 



(B8) 
(B9) 



In addition, one can show that ^34 = (^14. 

Notice that the combinations of helicity amplitudes contributing to the SSA, Eq. (I17p . are simply given by: 

|M ++;++ | 2 ee |i\7 1 T = 64. 2 a 2 e 2 ^ 



\M+-, + -\ 



|M 9 u | 2 = 647r 2 aV 



/ 2 



M+^+M*^^ = 64 7r 2 a 2 e 2 



s(-u) 



9 f2 



(BIO) 

(Bll) 
(B12) 



In the (transversity) g) (Collins) contribution to the SSA, the phase dependence of the last term above (<fii — = 
(f>— </>') combines with the k± phase in the transversity distribution (</>) and the Collins function phase (<pg), resulting 
in the simple expression given in Eq. (|17[) . 



Appendix C: Details for the computation of A^^ s 

In this Section we show some details of the explicit calculation of the transverse single spin asymmetry A^/^ s , 
Eq. (fl4")) . for the process £ — > hX, starting from the general expression for the polarized cross section given in 
Eq. (j4]). By performing the sum over all the helicity indices and taking into account that the helicity density matrix 
of a quark q can be written in terms of the quark polarization vector components, P q = (P£, Pjj, P?), as 

j/p.s (p\ + 4_ \ _\( 1+ p« ps \ (C1) 
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one obtains, for a spinlcss hadron h, 
E h & a {v,s)+i^h+x 

1 



d 3 P h 

' 2 {fq/P,s(x,k x ) (\M- 



16 7r 2 a; z 2 s 



d 2 ki d 3 p ± (5(p_L • p' ) J (p±) 5(s + i + u) 



(C2) 



++;++! 



|M_ 



+ ;-+ 



P« / g / p ,s(a:, ki) [Re(M ++;++ M* + ._ + ) cos$ - Im(M ++;++ M* + ._ + ) sin^] 
- P, 9 / g / P ,s(a:, [Im(M ++;++ M* + ._ + ) cos^ + Re(M ++;++ M* + ._ + ) sin^] 



(z,z»J-) } • 



In the above expression we have already extracted from the fragmentation functions D\ y (z,p±) their azimuthal 
dependence and exploited their parity properties (see Ref. [§0| for details): 



D++{z,Pi) = D—(z,p x ) = D h/q (z,pi) 
t) + _(z,p A _)=D + _(z,pi)e i < = ^A N D h/q 4z, P± )e< 
b + _(z,p i _) = [D_ + (z,p i _)Y . 



(C3) 
(C4) 
(C5) 



When computing the azimuthal asymmetry one has the difference of cross sections with opposite transverse spin, 
d<r(</>s) — d<r(</>5 + ir); using Eq. (flU]) and the definitions [8(| 



(C6) 



one obtains 



da(<f> s )-da(<f> s +w) = W d 2 kid 3 p ± S(p ± -p' q )J(pi)S( 



s + i+u) 



-Af q/pST (x,ki) (\M ++ ; ++ \ 2 + \M- +; -+\ 2 )D h/q (z,pi) 



A-f Sy/ST (x,ki)[MM ++ - ++ Ml + .,. + ) eos<-Im(M ++i++ M* + ._ + ) sin^] 



Af Sx/ST (x,ki)[lm(M ++;++ Ml + ._ + ) cos <f> q + Re(M_| 



) sin^]l A N D h/qt (z,pi) . (C7) 



Finally, using Eqs. ©, ©, flU]), (ESJ), ([BlOf - plli]) and the relations [8C 



A / Ss/Sl ,(a;,fe±) = 
A/ Sx/ST (a;,fc±) = 



h\(x, ki) 



hi(x, ki) + 



k 2 

^h 1T (x,ki) 
k\ L 



sm(<ps 



COS(0 S 



(C8) 
(C9) 



yields: 



da(^s) - d(j(0s + tt) = 51 / dx 2 dz 2 d 2 fc ± d 3 p ±( 5(p ± -p^) J(p ± )<S(5 + t + fi) 



A«f q/p i(x, ki) sin(0 s - 0) + |M 2 T) D h/q (z,p±) 



+ /» lg (x, fcj.) sin(0 s - 0' - <) MiM 2 ° A N D h/qt (z,pi) 



2M 2 



±S h±£(x, ki) sin(0 s - 20 + 0' + MfM" A N D h/qt (z,pi) 



(CIO) 



The first term on the r.h.s. of Eq. (|C10|) gives the Sivers contribution while the second term gives the transversity g) 
Collins effect. We have numerically checked that the third term gives negligible contributions. Notice that the various 
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terms of the type sin((/>s — $) appearing in Eq. (|C10|) can be decomposed as sin (j>s cos $ — cos</>s sin$: similarly to 
what has been explicitly shown in appendix I A 21 [see the comment after Eq. (|A27|) ]. the cos <fis terms integrate to zero. 
Thus, one obtains the simple expression of Eq. (fT7|l . given for (f>s = n/2. 
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